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Spinning compact binaries are shown to be chaotic in the Post-Newtonian expansion of the two 
body system. Chaos by definition is the extreme sensitivity to initial conditions and a consequent 
inabihty to predict the outcome of the evolution. As a result, the spinning pair will have unpre- 
dictable gravitational waveforms during coalescence. This poses a challenge to future gravity wave 
observatories which rely on a match between the data and a theoretical template. 
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Coalescing binaries are the primary objects of at- 
tention for future ground based gravity wave detectors 
such as LIGO and VIRGO. The successful detection of 
the waveforms requires a technique of matched filtering 
whereby the data is convolved with a theoretical tem- 
plate. Excellent agreement is required if a signal is to 
be drawn out of the noise. A possible obstacle to the 
method of matched filtering can surface if the orbits be- 
come chaotic. As shown here, the final coalescence of 
spinning, compact binaries proceeds chaotically for some 
spin configurations. Chaotic binaries with similar initial 
conditions may produce disparate waveforms and con- 
sequently they may not be detectable by the method of 
matched filtering. An alternative method must be sought 
for their detection. 

Many authors have emphasized that black holes are 
susceptible to chaos [0-|6|. Chaos has not received 
the attention it deserves in part because the systems 
studied have been highly idealized. An elegant ex- 
ample of chaos around black holes is provided by the 
Majumdar-Papapetrou spacetimes yj^l which arrange 
extremal black holes such that the gravitational attrac- 
tion of their masses is exactly countered by the electro- 
static repulsion of their charges. The spacetime is static 
and yields a simple solution. The geodesies however are 
formally non-integrable and fully chaotic Mm- A static 
spacetime produces no gravitational waves and so the 
chaotic scattering in the Majumdar-Papapetrou space- 
time remains just an interesting formal system, although 
gravity waves are produced by a third orbiting body [g . 
Chaos around Schwarzschild black holes has also been 
studied formally with a hypothetical perturbation of a 
test companion along the homoclinic orbits which mark 
the boundary between dynamical stability and instability 
M . Another important example of chaos around a black 
hole is the motion of a spinning test particle Q. This 
already shows the key features of the two-body system 
investigated here. 

In this paper, the most realistic description currently 
available of a black hole plus a companion is shown 
to succumb to chaos when the pairs spin. The Post- 
Newtonian (PN) expansion of the relativistic two-body 
problem P~[l2| provides the dynamical equations of mo- 
tion to 2PN-order [p3|jl^]. In the absence of spins, the 



existence of a conserved angular momentum and energy 
pot ensure that the swtem is in principle integrable to at 
least 5/2PN-order ||l^. The non-spinning pair still has 
two identifiable circular orbits for a given angular mo- 
mentum, one stable and one unstable. In the transition 
to chaos, the periodic orbits proliferate and these form 
the structure of the chaotic dynamics. The homoclinic 
orbits found in Ref. (l5| demarcate the region of phase 
space at which this occurs, perhaps at higher orders in 
the PN expansion. 

When spins are introduced at 2PN-order, the orbital 
plane precesses chaotically. There are now an infinite 
number of periodic orbits which form a fractal in the dy- 
namical phase space. We can isolate this fractal through 
the method of fractal basin boundaries H-0JI^-|I3| . Frac- 
tals are a particularly important tool in relativity since 
they do not depend on the coordinate system used, a 
point emphasized in [ [l8[ . 

In the notation of Ref. ||l^, the center of mass equa- 
tions of motion in harmonic coordinates are 



X = apN + aso + a-ss + apR- 



(1.1) 



The right hand side is the sum of the contributions to 
the relative acceleration from the PN expansion, from the 
spin-orbit (SO) and spin-spin (SS) coupling and from the 
radiative reaction (RR). The spins also precess by 

Si = (\iy.Si , 52 = fl2 X 5^2. (1.2) 

For brevity we do not rewrite the explicit forms of a 
and a here but they can be found in Ref. [Q. There 
are 12 degrees of freedom (x, a', 5i, 6*2). The form of 
eqn. (^_j) indicates that the magnitudes of the individual 
spins are conserved. To 2PN-order there is also a con- 
served energy E and a conserved total angular momen- 
tum J = L + S where L is the orbital angular momentum 
and S = Si + S2- In all, there are 6 constants of motion 
reducing the phase space to 6 degrees of freedom, plenty 
to allow for chaotic motion. The condition that the or- 
bit be perfectly circular r — r — (where r = \x\) still 
leads to an underdetermined set of equations for which 
there are an infinite number of spin configurations. This 
is evidence for the proliferation of periodic orbits and in- 
dicates the pursuit of an innermost stable circular orbit 
ill is futile. 
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1. The pair has mass ratio mijmx = 



4000 6000 

1.4/10 and no 

spins. The initial conditions are Xi/m = 10, yi = 0.3 and 
Zi — 0. Time is measured in units of the total mass m. Top: 
A 3D view of the orbit. Lower Left: The smooth phase space 
curve in the {r,f) plane. Lower Right: The waveform h+. 



Figure H^ shows typical orbital motion in the abs ence 
of spins and with the dissipative (RR)-term in eqn. (1.1) 
temporarily turned off. There is no precession of the or- 
bital plane and no chaos. Although the orbit is confined 
to a plane, the perihelion precesses within the plane due 
to the relativistic corrections. The regularity of the mo- 
tion is confirmed by the phase space diagram in fig. |l| 
which shows the motion to be confined to a smooth line 
in the (r, f) plane. The waveforms for specific orbits are 
obtained to 3/2PN-order using the results of Ref. |l^] 
and neglecting tail contributions. For simplicity we show 



the -|— polarization waveform, h^ — h^x, with the Earth 
located above the z-axis. 

If the compact objects spin, then the motion can be- 
come chaotic. The spin vector 5*1 is tilted by an angle 
9i measured from the z-axis and the spin vector 52 is 
tilted by an angle 02- The motion is clearly occupying 
three dimensions and is no longer confined to a plane 
as demonstrated in fig. g. A Poincare surface of section 
is constructed by plotting a point as the orbit crosses 
the z = plane from z>Otoz<0. A regular orbit 
would draw a smooth curve in the plane while a chaotic 
orbit speckles the plane with points unpredictably. The 
chaotic precession is indicated in the surface of section 
which has begun to turn to dust. The more titled the 
spin vectors, the thicker the dusty region in the surface 
of section. (Due to the large dimensionality of the phase 



space, the diagram is a projection onto the (r, f) plane. 
Cautious of any ambiguity this may introduce, we take 
the speckled surface only as confirmation of chaos seen in 
the precessional motion and the fractal basin boundaries 
discussed below.) The waveform is also shown. 

The binary of figure 2 could be a maximally spinning 
IOMq black hole with a rapidly rotating IAMq neutron 
star companion. The spins are each displaced from the 
initial orbital angular momentum by 45°. Large spin 
misalignments occur naturally in the formation of close 
black hole/neutron star pairs IpO]. The orbit shown is 
within the LIGO bandwidth with a frequency of roughly 
0(10 — 10^) Hz. With dissipation included, an orbit 
which begins regular at larger radii chaotically scatters 
as the pair draws closer and the signal sweeps through 
the LIGO bandwidth. 
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The initial conditions are 
Xi/m = 10, j/i = 0.3 and Zi = 0. The initial angles are 
9i ^ 02 = 45°. Top: A 3D view of the orbit. Lower Left: 
The surface of section in the {r,r) plane. Lower Right: The 
waveform h+. 

Chaos is not isolated to this specific binary. Instead 
of investigating individual orbits, we can broadly scan 
the phase space for chaos. There may be a sensitivity 
to the variation of any of the degrees of freedom as well 
as the relative masses of the compact objects. Since it 
is impossible to cover all variations, in this instance we 
limit our scan to search for chaos as the spin angles are 
varied. To do this, we look at a slice through the phase 
space which varies only the initial angle 9i of Si and the 



initial angle 62 of S2 for pairs which are otherwise given 
identical initial conditions (in this case TO2/TO1 = 1/3 and 
Si/m\ = 5'2/ml = 0.6). These could be black hole pairs. 
While the spins are consistent with neutron star pairs 
also, they are at such a close separation that tidal effects 
for the extended objects would be significant. The initial 
location in the (0i,02) plane is color coded black if the 
pair coalesce, grey if the pair separate by r/m > 1000, 
and white if stable motion is attained with more than 
50 orbits. A few pairs which separate to r/m > 1000 
may still continue orbiting. Increasing the cutoff would 
reduce the grey basin. Also, pushing the stable orbit 
condition to more than 100 orbits tends to increase the 
size of the black basins slightly as more orbits have a 
chance to coalesce. If there were no chaos, the boundaries 
between colors would be smooth while fractal boundaries 
signal chaos. The fractal basin boundaries of fig. Q clearly 
show a mingling of possible outcomes as the angles are 
varied. The extreme sensitivity to initial conditions is 
exemplified in the blown up regions in the lower panels 
of fig. which show the repeated fractal structure. 
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FIG. 3. Top: The fractal basin boundaries for pairs with 
m2/mi — 1/3 and Si/ra\ — 82/1^% — 0.6. AH orbits begin 
with Xi/m — 5,yi = 0.45. The initial angles [61,62) are 
varied. The axes are labelled in radians. 200 x 200 orbits 
shown. The middle and bottom panels are details of the upper 
panel. 

Compact pairs with initial conditions drawn from near 
the fractal basin boundaries will result in unpredictable 
outcomes. They will have correspondingly unpredictable 
waveforms. The waveforms for pairs selected from the 
initial conditions in fig. || are shown in fig. |4[ The orbits 
begin with nearly identical initial conditions. Although 
the difference in initial angles is only 3°, the waveforms 
are entirely different. The first pair separates while the 
second pair executes many thousands of orbits. 
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FIG. 4. The waveform h+ for pairs selected from the initial 
conditions in fig. H. Both orbits begin with d\ — 10°. The 
left panel began with 62 = 128° while the right panel began 
with 62 = 131°. The extreme angles were randomly chosen 
from the fractal set for illustration. Chaos is seen with more 
temperate angles as in fig. ^ 

It should be emphasized that orbits within smooth 
basins can still be chaotic. Well within the white stable 
basins, many orbits will precess chaotically as does the 
orbit of fig. 0. Similarly, many of the escape orbits and 
the merger orbits will chaotically scatter before reach- 
ing their final outcome. Fractal basin boundaries are a 
fairly blunt tool, insensitive to some manifestations of 
chaos. Therefore while fractal basin boundaries do prove 
the dynamics is chaotic, smooth basins are inconclusive. 

With the radiative reaction included, the pair goes 
from an energy conserving scattering system to a dis- 
sipative one. In any stability analysis, dissipation must 
be turned off to distinguish instability to the onset of 
chaos from instability to merger from simple energy loss. 
Once the chaos has been identified, radiative back reac- 
tion can readily be incorporated and we do so now. Under 
the effects of dissipation, some orbits will sweep through 
the chaotic region of phase space as they inspiral. The 
surface of section is not useful for a dissipative system 
since the radius of the orbit is shrinking as energy is lost 
to gravity waves. However, fractal basin boundaries are 
still effective at identifying extreme sensitivity to initial 
conditions. Another advantage is that several thousand 



orbits can be scanned at once. We use this method to 
show that dissipation does not obhterate the chaos. 

As energy is lost the binary pairs tend to coalesce in 
such a way that r — >merger is an attractor in phase space 
that can be described by another fractal set. To show 
this, we again look at an initial condition slice through 
phase space. We evolve each of these pairs under the 
influence of the radiative reaction force. We need to 
color code the initial conditions on the basis of some 
well defined outcome. Since all pairs considered coalesce, 
we have to select some other criterion than that used 
above. We choose to color code the initial location in 
the (^1, 62) plane white if the pair approach merger from 
below the z-axis and black if they approach merger from 
above the z-axis. The resultant fractal is shown in fig. 
0. Another criterion could have been selected and in this 
sense the basin boundaries are crude, as already men- 
tioned, but they are nonetheless powerful at signaling 
the presence of chaos. The conclusion to draw from this 
figure is that there is extreme sensitivity to initial spin 
angles for rapidly spinning, inspiralling IOMq black hole 
and 1.4Mq neutron star binaries. The pairs will inspiral 
along different paths as a result of this sensitivity and 
therefore will have disparate waveforms. Similar chaotic 
sets have also been found for different binary mass ratios 
and orbital parameters. 



would deter detection. Given that eccentricity in an oth- 
erwise simple orbit can greatly diminish the signal when 
matched against a circular template |gl[ , the chaotic pre- 
cession does not bode well. Still, the luminosity in grav- 
ity waves is enhanced for some of these wilder orbits m , 
as was already seen along the regular homoclinic orbits 
p5[ . Though unlikely, an optimist might hope that di- 
rect detection of these gravity waves will be possible if the 
signal is boosted substantially above the noise, relieving 
the dependence on a theoretical template. 

The inherent difficulty in the direct detection of grav- 
ity waves highlights the importance of indirect methods 
of detection. Corroborating evidence for gravity waves in 
electromagnetic observations may be promising. Chaos 
can have unexpected benefits if the black hole is able to 
capture the light from a luminous companion for many 
chaotic orbits before some of the light escapes. Such 
chaotic scattering of a pulsar beam around a central black 
hole could lead to a diffuse glow around the pair [^. 
While this signature is likely to be faint, any confirmation 
of a gravity wave signal will be welcome. 

I am grateful to E.J.Copeland, R. O'Reilly, and 
N.J.Cornish for their valuable input. This work is sup- 
ported by a PPARC Advanced Fellowship. 
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5. The fractal basin boundaries with dissipation in- 
The parameters are 7712/7711 = 1.4/10 and spins 



nil a-nd S2 ~ 0.7 m^. The orbits begin with 



FIG. 
eluded, 
Si = 

Xi/m — 26, yi — 0.15 and Zi — 0. The pair can execute 
anywhere from to 0(40) orbits before coalescence. The ini- 
tial angles (61,^2) are varied from —n to n. 300 x 300 orbits 
shown. 

This work demonstrates the existence of chaotic re- 
gions of phase space. At least some orbits will move 
into this chaotic region as they inspiral. Of course some 
orbits will still be regular such as circular inspiral with 
spins exactly aligned with the orbital angular momen- 
tum. A systematic scan of all parameters is needed to 
ascertain when the dynamics is predictable and regu- 
lar and when it is chaotic. A quantitative comparison 
of the waveforms from a chaotic orbit against a circular 
template is also needed to evaluate how seriously chaos 
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